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Heavy-ion collisions at the BNL-RHIC collider can probe 
whether gluon saturation effects in nuclei at small x have set 
in, or whether leading-twist perturbative estimates of particle 
production are still applicable. I discuss that soon to come 
data from pA collisions at RHIC may provide more systematic 
insight into this problem than present data on the energy and 
centrality dependence of hadron production from Au + Au 
collisions. Results from pA collisions at RHIC could also help 
to understand some controversial and puzzling results from 
ion-ion collisions, for example the large azimuthal asymmetry 
at high p±, as measured by the STAR collaboration. 



I. INTRODUCTION 

Perturbative (short-distance) QCD predicts a double- 
logarithmic distribution of gluons in a hadron or nucleus 
at small x (i.e., at rapidities far away from the valence 
quarks): dN g ~ dk±/k~ L dx/x. This result, however, is 
derived in the dilute limit where some effects, like gluon 
recombination, have been neglected. Therefore, as k± 
decreases, the DGLAP approximation should eventually 
become invalid, and the predicted rise of the gluon den- 
sity should flatten out. 

Namely, at large distances (though still much smaller 
than the confinement scale) the non-linear terms in the 
Yang-Mills equations tame the growth of the color charge 
density in hadrons as compared to the power-law 
growth predicted by the linearized DGLAP theory. The 
gluon density should saturate at distances where the color 
charge density p a per unit transverse area becomes of 
order l/<? 2 , where g is the coupling, and where there- 
fore QCD-evolution becomes genuinely non-perturbative. 
Such large occupation numbers are characteristic of a 
condensate, and so classical methods may be useful in de- 
termining the color field in a hadron or nucleus at small 
x, i.e. at high density (the "saturation effect" is 

expected to show up in large nuclei much "earlier" , at 
larger x, than for protons). The local color charge den- 
sity is a stochastic variable (fluctuating randomly as one 
moves around in the transverse plane) and observables 
eventually have to be averaged over some distribution 
of p a . This is analogous to spin glasses, and hence the 
name "Color Glass Condensate" g for the structure of 
a nucleus at small x. 

The Relativistic Heavy-Ion Collider (RHIC) at BNL is 
dedicated to the study of high-energy hadronic interac- 
tions, from p + p over p + A to A + B collisions. One 



of the important issues to settle is whether the above- 
mentioned gluon saturation scale Q s at the values of x 
that are relevant at RHIC energy is in fact (substantially) 
larger than the QCD confinement scale A. Optimistic es- 
timates || do indeed yield Q 2 s (y = 0) ~ 2 GeV 2 for 
central Au+Au collisions at RHIC, that is Ql/A 2 ~ 50. 
If so, then perturbative estimates of particle production 
at leading twist Q can not be employed to compute the 
dominant part of the inelastic cross section because one 
needs to resum all powers of a s times the charge density 
squared per unit of transverse area ( "higher twists" ) . In 
other words, particle production will be dominated by the 
decay of the classical non-abelian field, except at very 
short distances (truly high p t ) where gluon occupation 
numbers are small and one enters the region of applicabil- 
ity of perturbation theory. The initial conditions for the 
subsequent evolution, which possibly proceeds through 
a so-called Quark-Gluon Plasma, then have to be ob- 
tained from the Color Glass Condensate model ||. On 
the other hand, if it turns out that at RHIC energy Q s 
is not substantially larger than A, as some conservative 
estimates indicate |J , then this means that leading- twist 
perturbation theory is safe at RHIC, and that it can be 
employed to describe the bulk of particle production . 
To understand particle multiplicities and energy densi- 
ties in the central region one does not need to consider 
the soft classical field then Q. 

Naturally, as RHIC started operation with collisions 
of gold ions, theoretical studies focused on the energy 
and centrality dependence of particle production from 
the CGC model versus that from leading twist perturba- 
tion theory with fixed infrared cutoff p Q ~ 2 GeV ||,§,|l0| ■ 
Within the CGC model, a rather mild increase of 
dN c h/drj both with centrality (or the number of partic- 
ipants) and with energy is expected. In contrast, phase 
space for particle production at p± > po rapidly opens 
up in perturbation theory, and so the growth is faster. 
The energy dependence of particle production at RHIC 
turns out to be quite weak pi| , approximately consistent 
with saturation models |p|,[l0[ However, to my opinion 
it is hard for both saturation models and soft-hard two- 
component models to make absolute predictions on the 
10-20% level (at present), which would be required to fal- 
sify one or the other model using just the limited energy 
range of RHIC. Similar concerns apply to the centrality 
dependence of particle production [Q . This observable 
is plagued in particular by the uncertainties in determin- 
ing experimentally the number of participants for large 
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impact parameter collisions. In that respect, it would be 
cleaner to determine dN c h/dr]/N part by using only cen- 
tral collisions but varying the mass number of the collid- 
ing nuclei. Of course, this will not be available within the 
near future, although it is doable in principle at RHIC, 
which can collide essentially any ions from p to Au. 

Is it possible then to verify or falsify the gluon sat- 
uration in nuclei at RHIC energies ? I think a good 
opportunity for that exists in pA since protons are likely 
quite dilute at RHIC energy, while saturation may occur 
for large nuclei |l3|,|l4j]. In that case, distinct signatures 
emerge in the transverse momentum and rapidity distri- 
butions of secondaries, which are correlated (and so can 
be tested experimentally without necessarily relying on 
absolute predictions). Also, one might be able to test 
various scalings of observables like dN/dy, (p±) and so 
on. More detailed observables like the p± -broadening of 
minijets at forward rapidities fl5|| , or diffractive single-jet 
production |lq| may as well provide deeper insight. 



II. PROTON-NUCLEUS COLLISIONS 

pA collisions at high energy are very interesting be- 
cause one can actually compute the total multiplicity 
analytically |0,0|, and also make some predictions for 
transverse momentum spectra and the scaling of (p±) 
with the produced multiplicity dN/dy. 

The key point is that the color charge density per 
unit transverse area is much smaller in a proton than 
in a big nucleus. Therefore, the saturation momentum 
of the proton is much smaller than that of the nucleus, 
Q { s 1} < Q ( s ] - At p± > Qi 2) the fields of both the proton 
and of the nucleus are weak, and so one expects that or- 
dinary perturbation theory is applicable. In particular, 
to compute single-inclusive gluon production at high p± 
the well-known expressions from collinear factorization 
(with DGLAP evolution) can be employed. 

On the other hand, the inelastic cross section (particle 
production) is dominated by the kinematic region 

(2) 

p± < Qs ■ In that regime, the field of the proton is still 
weak, and so can be treated within perturbation theory; 
but the field of the nucleus is in the non-linear regime, 
and one must resum interactions of the radiated gluons 
with the background field of the nucleus to all orders. 
A simple way to understand why is to realize that while 
each additional gluon "line" from the produced gluon to 
the nuclear field comes with an additional power of the 
coupling a s , the occupation numbers of the nuclear field 
at (or below) the scale Q s are of order l/a s , and so 
additional "rescatterings" are not suppressed by powers 
of the coupling constant. 

The calculation of gluon production essentially consists 
then of solving the Yang-Mills equations in the forward 
light-cone to all orders in the background field of the nu- 



cleus, but to leading order in the field of the weak proton. 
From the Fourier transform of the field at asymptotic 
times one reads off the amplitudes of the modes, which 
are then squared and averaged over the color charge dis- 
tributions in the proton and in the nucleus, using a Gaus- 
sian weight |l| . Finally, this leads to the gluon distribu- 
tion El 
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the gluon propagator. Xiiu) denote the effective den- 
sities of color charge as "seen" by the gluon produced 
at rapidity y. It is the integrated charge density of ei- 
ther one of the sources, from its beam rapidity down to 
y. Note that (|l|) is linear in \i but resums all orders 
in xi- A numerical evaluation of the above result has 
not yet been performed, mainly because one needs to de- 
termine the dependence of \i on rapidity 13]. One can 
nevertheless gain some insight by considering two lim- 
its. Namely, at asymptotically large k±, the exponential 
can be expanded to first non-trivial order, leading to the 
well-known pQCD result 
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where we introduced the shorthand notation Q s ' (y) = 
^/g 4 iV c xi :2 (j/)/87r |T7| . Thus, one recovers the stan- 
dard perturbative ~ a^/k\ behavior at very high k±, 
with a logarithmic correction analogous to DGLAP evo- 
lution. Note that xi, X2 scale as A 1 / 3 and AtJ 3 Q, re- 
spectively iQ, while the integral over d 2 b gives a fac- 
tor of ttR 2 oc A^ . Therefore, in this kinematic region 
dN/d 2 k±dy scales like A\^ 3 A2, up to logarithmic cor- 
rections. This holds also for the integrated distribution 
dN(k± > po)/dy above some fixed ^-independent scale 

Po, as in the soft-hard two-component model. On the 

(2) 

other hand, when integrating over k± from Q s to infin- 
ity, the contribution from large k± to the rapidity density 
is 
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Again, the integral over d 2 b gives a factor kR\ ~ A^ 3 , 
and so dN/dy scales like A^A 2 , . 

On the other hand, when Qi X \y) <C k± <C Qs (y), 
one finds that 
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where £(y) = Qs (y)/Qs' > (y)- From dimensional con- 
siderations, it has been suggested |2(| that in symmetric 
A+A collisions, and at central rapidity, (k±) 2 scales with 
the multiplicity per unit of transverse area and of rapid- 
ity, (k±) 2 cx dN /d 2 bdy. A similar scaling relation can be 
derived from eqs. fi[7fi for the asymmetric case, 



This form ~ agXi/k^ is to be compared with that from 
ec l- (H)i ~ a lxiX2/k 4 L , valid at high k±. A schematic 
distribution in transverse momentum is shown in Fig. [l]. 
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Thus, (k±) 2 is proportional to the multiplicity per unit of 
rapidity and transverse area, times a function of the ratio 
of the saturation momenta. If source one is very much 
weaker than source two, i.e. in the limit | log£| ^S> 1 — £, 
the third factor on the right-hand-side of (@) depends on 
log£ only. Neglecting that dependence, and assuming as 

1/3 

before that xi,2 are proportional to A x 2 , one has the 
approximate scaling relation 
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FIG. 1. Schematic k± distribution for gluons produced 
in high-energy A\ + A 2 collisions at rapidity y such 
that Qs (y) "C Qs (y)- In the perturbative regime, 
dN/dk±dy ~ 1/^1- In between the saturation scales for the 
two sources, dN/dk±dy ~ 

Although the figure does not provide a quantitative 
prediction for experiments at RHIC, one should notice 
non-trivial features. It is predicted that when strong- 

(2) 

field effects set in around Qs , the k± distribution of 
secondaries flattens. Not only does the scaling with A 2 
differ in the two regimes, as mentioned above, but more 
importantly the "turnover" point Q s (y) is a function 
of rapidity ! Experimentally, one can thus take the ra- 
pidity dependence of the saturation momentum from a 
parametrization of HERA data which also seems 
to fit the observed rapidity dependence of dN/dy from 
Au + Au at RHIC || , and test whether the turnover in 
the transverse momentum distribution from fig. [l] moves 
as one changes the rapidity in a way consistent with those 
parametrizations . 

From (g), the fc^-integrated multiplicity in the nonper- 
turbative regime Q s (y) < k± < 
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Thus, at fixed impact parameter, the multiplicity scales 
as A\^ 3 , up to the square of a logarithm of (A2/A1) 1 / 3 . 

The average transverse momentum in the saturation 
regime is given by 

e-i-iogc 



(k±) = 2 Q 
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In practice though one expects significant corrections to 
the simple scaling relation (^), as given by eq. (|sj) - pA 
collisions provide a natural testing ground for the scal- 
ing of (p±) with dN/dy since distortions of the above 
scaling relations from final-state interactions ( "hydrody- 
namic expansion") should not be a major issue. 

TABLE I. Scaling of various quantities in the perturbative 
and saturation kinematic regimes. 



sat. 



pert. 



dA^/d 2 6d 2 fc ± d?/ A\ /z 
(k ± ) 2 (AiA a ) 
A. 1™ dN 

dy 1U & d'-'k ± dy 



A\ /3 Al /3 



~ 1/ ' 3 dN/dy 1 (for fixed po scale) 



d 

< I'.' 



logQ 3 1) (y) + \ogQ 3 '"(y) 



(2), 



3 



For rapidities far from the fragmentation region of 
the large nucleus, and for Qs(y) <C k± <C Qs (y), 
dN/d 2 k±dy varies with rapidity like 

d diV 



« 9 x'i(y) 
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dy d 2 k±dy 

Thus, an experimental measure for the evolution of the 
CGC density parameter in rapidity is 

dlogdV/d^dy dlogxi(y) 



dy 



dy 
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Now consider the case of high k± ^> Q s (y) described by 
eq. (|^). In that regime the rapidity distribution is pro- 
portional to Xi(y)X2(y), an d so dN/dyd 2 k± varies with 
rapidity like 



dlogd7V/d 2 fc±dy _ dlogxi(y) , dlogx 2 (y) 



(12) 



dy dy dy 

Subtracting (pd| ) from ( p"2|) provides an experimen- 
tal measure for the renormalization-group evolution of 
X2(y)- This does not only provide an independent verifi- 
cation of the HERA fits , but also a consistency check 
on how the turnover point (where the distribution flat- 
tens) from fig. [l| moves as one goes to different rapidity ! 
In table 1 1 have listed a few scaling relations of various 
quantities. 

In fact, the results obtained above may even be rel- 
evant for collisions of equal-size nuclei, like Au + Au. 
Namely, if we consider a collision at non-zero impact pa- 
rameter, as shown in fig. 0, it is easy to realize that near 
point A in that figure the situation is precisely the one 
considered above for pA collisions: a low-density projec- 
tile colliding with a high-density target. 




Target 



FIG. 2. Illustration of a collision of equal-size nuclei at 
non-zero impact parameter. Near point B, close to the edges 
of both nuclei, the local densities of partons in the initial state 
are small; there, the conditions resemble those for pp colli- 
sions. Near point A, we are encountering a pA-like situation: 
the local parton density near the edge of the projectile is small 
but it is large in the target (close to the center). 



From the discussion above, the k± distribution of glu- 
ons produced in the impact is harder near point A than 
near point B, due to saturation effects in the target. Now, 
if those gluons suffer from radiative energy loss as they 
interact with each other, it is clear that a detector posi- 
tioned at the top of the page will detect less momentum 
flowing into it than a detector positioned at the right of 
the page. Thus, radiative energy loss will generate an 
azimuthal asymmetry 
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at transverse momenta on the order of a few GeV [ED. 
The point here is that v% will be larger if the p±- 
distribution of produced gluons is harder at A than at 
B as compared to the case when diV/d 2 p_L is the same 
at A and B. This might be important for understand- 
ing the large values of V2 measured by STAR at RHIC at 
p± — 2 — 6 GeV [Q , as it has been pointed out that large 
energy loss by itself (assuming the same p^-distribution 
throughout the overlap zone in fig. |J) does not lead to 
large enough V2 p3f , in particular for semi-central events 
(10% centrality class). Quantitative computations of 
V2 from the gluon saturation model are presently under 
way || . 

So far, all attempts to understand the large v% at high- 
p± invoke radiative parton energy loss p5fl 



. 3a Srl , 2E 
AE ~ — -E cr log - — j 



(14) 



Here, L is the typical length of the medium traversed 
by the jet, mg is the IR cutoff provided by the mass of 
static electric gluons, E > E cr is the jet energy, and 
E cr = m 2 L 2 /X. AE/E depends on the density of the 
medium via the mean-free path A of the radiation. Going 
down in -^/s, presumably this should reduce V2 at high 
pj_ (in the sense that the leading-twist prediction V2 = 
should be approached faster). Is this the case ? 

Finally, if indeed V2 at p± = 2 — 6 GeV is generated by 
final-state interactions of the minijets, it might be sensi- 
tive to nearly critical scattering at T c p6| , which would 
evidently be very exciting. QCD with 3 colors is not 
an ideal gas near T c ; rather, the pressure p(T) is signif- 
icantly smaller than the ideal-gas pressure Pid(T). The 
transition from confinement to deconfinement appears to 
be nearly a second order transition [p7| , in which case 
(mi/T) 2 ~ y/pjpid drops rapidly as T c is approached 
from above ]26fl . This behavior is quite different from a 
strong first-order transition, and can enhance energy loss 
(and V2) near T c . 



III. CONCLUSIONS 

Au + Au collisions from RHIC may have offered a first 
dance at gluon saturation effects in large nuclei: exper- 
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imentally, the energy and centrality dependence of par- 
ticle production is quite weak. However, in my opin- 
ion it is too early yet to draw definite conclusions from 
the Au-runs at ^ = 130 AGeV and 200 AGeV. Col- 
lisions at lower energy and in particular with different 
projectile/target combinations should be analyzed in the 
future. Collisions of small projectiles on large targets, 
for which I use the generic acronym "pA" , could reveal 
distinct scaling laws predicted by the Color Glass Con- 
densate model. 
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